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The potential of constraints on the angle a of the Unitarity Triangle from studies of — > n + K~7t° decays is reviewed. The experi- 
mental inputs needed for an isospin analysis are starting to become available. The precision of current measurements is extrapolated to 
higher luminosities and the constraints on a that could be obtained in the foreseeable future are discussed. Studies of the sensitivity of 
a time-dependent Dalitz plot analysis are briefly reviewed. 



1 Introduction 

The BABAR and Belle Collaborations have performed 
searches for CP violation in Z?°/B° decays to % + %~ [[Q, 
where the mixing-induced CP asymmetry is related to the 
angle a of the Unitarity Triangle (UT). The implications 
of these results are discussed in the contribution [ [2) to 
the present proceedings. In this paper, we discuss the 
K + K~K° final state which probes both direct and mixing- 
induced CP violation. The mixing-induced CP violation is 
related to a. 

The decay Z?°/B° — > 7C + 7r~7r° is dominated by the res- 
onant intermediate states p + n~ and p~n + [ |3] |4j 0. 
The leading Feynman diagrams for the B° — > p ± n zf de- 
cays are shown in Fig. 0a). Note that Fig. [Ha) repre- 
sents two distinct diagrams, one where the p meson is 
formed from the W + boson, and one where the p meson 

is formed with the spectator quark. As pointed out in [|6j, 

o — 

these diagrams mean that both B and B can decay to the 
p + n~ state, and therefore that CP violation can manifest 
itself in the interference between decays with and without 
mixing. The moduli of the two corresponding amplitudes 
are expected to be comparable in size, and CP-violating 
effects could thus be large in these decays. The same com- 
ment applies to the p~n + state. Contributions from pen- 
guin diagrams (see Fig.^b)) could be sizeable and com- 
plicate the interpretation of measured CP asymmetries in 
terms of the elements of the CKM matrix. The decay 
B°/W — > % + % n° can also proceed through the p°n° in- 
termediate state. The corresponding (colour suppressed) 
Feynman diagram is shown in Fig. ^c). This decay has 
not been observed yet; current experimental limits are of 
the order of 5 • 10~ 6 [|3]|4][5). Contributions from scalar 
resonances like, e.g., /o(980) — > and higher exci- 

tations are also expected. Non-resonant contributions to 
B°/W — > K + K 7T° are small, they are found to represent 
less than 4 % of the total rate [[5). 

Different approaches to the analysis of — > n+n^n® 
are used in the literature. In the so-called quasi-two-body 
approach [ 7 1 one restricts oneself to the two regions of the 




d d 

Figure 1. Feynman diagrams fori? — * p ± 7t^ : andB — * p°7l°. 



n + n~n Dalitz plot dominated by either p + n~ or p~n + . 
Another approach is the time-dependent Dalitz plot anal- 
ysis originally proposed in [[5]. It is optimal in the sense 
that it exploits all available information, but experimentally 
much more involved. Specifically, the full dynamics of the 
K + K~K° final state are studied, and the interference be- 
tween the p + 7t~, p~7t + and p 7T° contributions in differ- 
ent points of the Dalitz plot is used to obtain simultaneous 
constraints on the weak and strong phases in these transi- 
tions. This analysis allows the measurement of a even in 
the presence of non-negligible penguin contributions. 

The quasi-two-body approach is discussed in section |2] 
Recent BABAR results obtained in that framework are re- 
viewed. Based on this and other experimental results, the 
sensitivity to a is extrapolated to higher luminosities ex- 
pected to be available in the future. Present data samples 
are not large enough for a time-dependent Dalitz plot anal- 
ysis. The expected statistical sensitivity and some of the 
major difficulties expected in this analysis are briefly dis- 
cussed in section|3] 
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2 Quasi-two-body analysis 

The observables in this analysis are defined in section [2~TI 
A recent BABAR measurement of these observables is sum- 
marised in section l2~2l In section l2~3l the present sensitiv- 
ity to a is assessed and extrapolated to higher luminosities. 
The discussion in sections |2~TI and 1231 follows essentially 
the discussion in [ 9 1, and the reader is referred to this paper 
for additional details. 

2.1 Basic formulae and definitions 

The four decay amplitudes relevant for the quasi-two-body 
analysis can be written as 



(1) 



where the T and P symbols represent the effective com- 
plex tree and penguin contributions. Contributions with 
other weak phases are neglected here. The parameters R u = 
\VudV: b \, R c = \V cd V* b \ and R t = \V td V*\ are the sides of 
the UT, /3 and y follow the usual convention for the angles 
of the UT The unitarity relation R u e iy + R c e in + R,e'P = 0, 
valid to €?(X 5 ) in the Wolfenstein expansion parameter A, 
has been used to express the contributions from c quark 
loops in the penguin amplitudes in terms of the other con- 
tributions. 
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where the parameters p and q describe the composition 
of the B° mass eigenstates in terms of the flavour eigen- 
states [ 7|. These definitions are inspired by the definition 
of the parameter A that is commonly used in the discus- 
sion of B decays to CP eigenstates (see, e.g. [ 7|). Each of 
the two A H (~ + ' parameters defined above involves only 
one pn charge combination. They are insensitive to di- 
rect CP violation, but their imaginary part is the sum of 
contributions from the weak phase a and strong phases. 
The ff + ~' h ) parameters involve both pn charge combi- 
nations, but only one of the two amplitudes represented 
in Fig.^a), i.e. either corresponding to transitions where 
the p is formed from the W or from the spectator quark. 
Their moduli are linked to direct CP violation, while then- 
phases measure effective weak angles oc^ ' + ' defined be- 
low. The parameters A + ~(~ + ) do not have the same prop- 
erties as the parameter X in the analysis of CP eigenstates, 
i.e. A + ~ 7^ ±1 or X h ^ ±1 does not necessarily imply 



CP violation. We define the parameters Acp and A tag such 
that 
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Acp has the desired properties, i.e. Acp ^ 1 implies CP vi- 
olation. Atag is a measure of the "self-taggedness" of the 
B Q /B° -> p^n* decay. If, e.g., A+~ =A~ + = 0, i.e. the 
charge of the p identifies the flavour of the B (the B — > p n 
decay is self-tagging), then A tag = and mixing-induced 
CP violation cannot occur. For the other extreme case, 



\A~+\ and |A" 



\A + |, we obtain |A tag | = 1. 



At an asymmetric B-factory like PEP-IIABABAR, the time- 
decays is studied using 
Z?°S°, where one of the 



dependence of B° /B° 



P^TZ 



pairs of B mesons from T(4S) 
B mesons (denoted B pn ) is reconstructed in the pn final 
state, and the flavour of the other B meson (denoted Stag) 
is determined (tagged) from its decay products. Defining 
At = t„ k ~ f tag as the proper time interval between the de- 
cays of the B pjz and the B ta g, the time-dependent decay 
rates are given by 
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1 + QtegiSpx ± AS pK ) sin(Am d Af ) 
- Gtag (C p n ±AC pn ) cos (Am d At) 

where gtag = 1(— 1) when the tagging meson B® ag is a 
B°(B°), X is the mean B° lifetime, and Am d the mixing 
frequency due to the eigenstate mass difference. The ob- 
servables Aqp, S pK , ASpx, 



terms of the amplitudes defined in Eq. Q: 



Cp K and ACpj: can be written in 
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The quantity A^p measures time-integrated and B-flavour- 
integrated direct CP violation. Flavour-dependent CP vi- 
olation is parameterised by C pn . S p7C measures mixing- 
induced CP violation related to the angle a. The quantities 
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ACpn and AS pn are insensitive to CP violation. AC p7C is 
a measure of the self-taggedness of the B° /B° — > p ± n* 
decay, and AS pjl is related to the strong phase differences 
between the different T and P amplitudes in Eq. Q. 

We adopt the definition of the effective weak angles 
Cttff and a eff + given in Ref. r fTol : 

2a e +" = arg jc~ + , 2a eff + ee arg jc+~ . (5) 

The relation between the observables in Eq. and these 
effective weak angles can be written as 

S pK + ASpn = Tjl-iCpa + ACpx) 2 sin(2a e +- + 5) , (6) 
S P n - AS pn = ^-{Cpn-ACpx) 2 sin(2a eff + - 8) , (7) 

where 8 = arg(A h A H *). The parameter 8 is related 
to the phase between the amplitudes A(B° — > p + Tt~) and 
A(B° — > p~7£ + ). It cannot be determined in a quasi-two- 
body analysis as one restricts oneself to regions of the 
Dalitz plot that are dominated by either p + n~ or p~n + . 
Information from the region of the Dalitz plot where p + n~ 
and p~% + interfere is needed to determine 8. In the case 
of negligible penguin amplitudes, the effective weak angles 
are equal to the angle a of the UT. 

2.2 BABAR measurement of CP- violating asymmetries 

in B° -> p ± n^ and B° -> p K+ 

A first preliminary measurement of the parameters A^p, 
Sp K , AS px, C p7[ and AC pK has been presented by the BABAR 
Collaboration at ICHEP 2002 [ II II . An improved version 
of the analysis with significantly reduced systematic uncer- 
tainties has recently been finalised N12I . This update also 
includes precise measurements of the branching fractions 
= ,@{B°/B° -> p ± n*) and 39$ = ,%{B°/B Q -> 
p±^). 

BABAR' s internally reflecting ring-imaging Cherenkov de- 
tector provides good separation between kaons and pions 
from B — ► p% and B — > pK (the typical separation de- 
pends on the particle's momentum and varies from 8(7 
at 2 GeV/c to 2.5a at 4 GeV/c). To treat the small 
overlap at higher momenta in an optimal way, a simul- 
taneous analysis of the p% and pK states is performed. 
The time-dependence of the B — > pK decays is given by 
Eq. (0 with CpK = S p k = AS p k = and AC p k = — 1 as 
the B — > pK decay is self-tagging. As the pn/pK branch- 
ing fractions are small [£?(10~ 5 )] and the p resonance is 
broad, the pit and pK channels suffer from large back- 
grounds. Continuum e + e~ — > qq(q = u,d,s,c) events are 
the dominant source of background. As these events are 
very abundant, their properties can be extracted from the 
data simultaneously with the parameters of the pn/pK sig- 
nals. Cross-feed from other B decays (both charmless 



and b — > c decays) is potentially more dangerous, as some 
of the corresponding final states are very similar to the 
signal, and these backgrounds are hidden, together with 
the signal, in a much larger number of continuum back- 
ground events. B candidates are identified using the beam- 
energy substituted mass mgs [ 1 1 1 and the difference AE 
between the centre-of-mass (CM) energy of the B candi- 
date and where is the total CM energy. AE also pro- 
vides good discrimination between signal and background 
from higher-multiplicity B decays, whose branching frac- 
tion is not well known in many cases. To enhance the dis- 
crimination between signal and continuum, a neural net- 
work (NN) is used to combine four discriminating vari- 
ables: two event-shape variables (continuum events tend 
to be jet-like, while B events are more spherical), the mass 
of the p candidate, and cos 9 K , where 9 K is defined as the 
angle between the n° momentum and the negative B mo- 
mentum in the p rest frame. The pit and pK event yields as 
well as the observables discussed above are determined in a 
simultaneous likelihood fit. The event variables used as in- 
put to the fit are mgs, AE, NN output, n/K Cherenkov an- 
gle and At with its uncertainty a At- The likelihood function 
contains contributions for signal, continuum and B-related 
backgrounds. 24 parameters that provide an empirical de- 
scription of the continuum events are also free to vary in the 
fit. Cross-feed from other B decays is studied using simu- 
lated events. The branching fractions of unmeasured decay 
channels are estimated within conservative error ranges. 
More than 100 charmless decay modes have been studied, 
and the 3 1 which contribute to the final event sample, plus 
the contribution from b — > c decays, are taken into account 
in the likelihood function. Backgrounds from two-, three-, 
and four-body decays to p % are dominated by B + — > n + n°, 
B + — > p°n + , and longitudinally polarized B° — > p + p _ 
decays, respectively. The pK sample receives dominant 
two-body background from B + — * K + n° and three- and 
four-body background from B — > K*n and higher kaonic 
resonances, estimated from inclusive B — > Kkk measure- 
ments. Using 89 x 10 6 BB pairs collected at the T(45), 
the BABAR Collaboration finds 428 ± 34 (120 ± 21) p ± n^ 
(p ± K^ : ) events and measures H12I : 

SSp = (22.6 ±1.8 ±2.2) x lO" 6 , 
= (7.3 + ii±1.3)xl0- 6 , 

Ag = -0.18 ± 0.08 ± 0.03 , A£p = 0.28 ± 0.17 ± 0.08 , 
C p% = 0.36±0.18±0.04, S p% = 0.19 ±0.24 ±0.03, 
AC pn = 0.28+o;|g±0.04, AS p7I = 0.15 ±0.25 ±0.03 , 

where the first errors are statistical and the second sys- 
tematic. The dominant systematic uncertainties are due 
to the uncertainties in the model of the B-related back- 
grounds, and possible interference between the doubly- 
Cabibbo-suppressed b — > Tied amplitude with the Cabibbo- 
favoured b — > cud amplitude for tag-side B decays [ 1131 . 
Other important uncertainties, especially in the branching 
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fraction measurements, are due to the uncertainties on the 
shapes of the reference distributions for signal used in the 
likelihood function, and to fit biases caused by small ef- 
fects that are not included in the likelihood model. While 
the branching fraction measurements are limited by sys- 
tematic uncertainties, the other measurements will remain 
statistics limited for quite some time. The uncertainties 
due to B-related backgrounds are expected to improve as 
more experimental constraints on the corresponding decay 
modes become available. These results are consistent with 
CP conservation at the 2.3(7 level. 

2.3 Expected constraints on a 

Figure |2ja) shows the constraints on a obtained using the 
BABAR measurement and the assumption that penguin am- 
plitudes are negligible (P + ~ = P h = 0). The eightfold 
ambiguity is due to the presence of the unknown strong 
phase 8. This figure also shows the constraints obtained 
when, in addition, the strong phase 8 is assumed to be zero 
(P = 8j — 0). This exercise is purely academic as the pen- 
guin contributions are not expected to be negligible and the 
strong phase 8 = arg(T h T + ~*) is unknown. It is shown 
here to illustrate the accuracy that could be obtained if the 
penguin amplitudes and the strong phases of the tree dia- 
grams were known (e.g. from theory) and to set the scale 
for the discussion that follows. In this case, a could be 
determined up to a twofold ambiguity with an accuracy of 
7° per solution. This is better than the indirect constraints 
on a from the standard global fit of the CKM matrix (see 
Fig.Ea)) [M- 

With the assumption of SU(2) flavour symmetry, the pen- 
guins can be constrained using experimental information 
on the other B — ► pn charge combinations. The key ob- 
servation is that gluonic b — > d penguins are pure AI = A, 
while tree-level b — > uud decays have AI = | and AI = j 
components. With the definitions A +0 = A(B + — ► p + 7T°), 
A°+ = A(B+ -> p°7T+) andA 00 = A(B° -> p°7Z°) this leads 
to the two pentagon relations 

V2(A+°+A 0+ ) = 2A 00 +A+-+A-+, (8) 
V2{A +0 +A 0+ ) = 2A° +A + -+A-+ ; (9) 

see [ 7j for the details of the isospin decomposition. The 
effects of electroweak penguins and electromagnetic tran- 
sitions, which can have AI — f components, are neglected 
here. This leads to 12 unknowns (6 complex amplitudes 
T + ~, P+-, T-+, P +, T 00 and T +0 , and a, minus one ar- 
bitrary global phase) and 13 observables (see Table^. The 
isospin analysis is not yet feasible with present statics from 
the B-factories (some of the relevant decay modes have not 
even been observed, Table [Q. A rough extrapolation into 
the future is presented in [ 9 1 . With an integrated luminos- 
ity of 500 fb~ 1 the experimental situation might look like in 
Table|2] One of the infinite number of sets of values for the 



12 unknowns that reproduce the BABAR measurement has 
been chosen to obtain the central values. The uncertainties 
are educated guesses based on the uncertainties of present 
measurements and searches, (optimistically) assumed to 
scale with 1 / \/N. The constraint on a obtained in this 
scenario is shown in Fig.[2b). This figure also shows the 
constraint for even higher luminosities of 2 ab~' (achiev- 
able by the first generation B-factories ?) and 10 ab _1 (Su- 
per B-factory). Very high luminosities are needed to obtain 
meaningful constraints on a, and even at high luminosi- 
ties several solutions exist. Figure |2jc) shows a slightly 
different scenario where ,<%(B° — > p 7T°) is assumed to be 
significantly smaller than in the first scenario (below the 
experimental sensitivity). In this case, the SU(2) analysis 
provides meaningful constraints above 2 ab~'. 

3 Time-dependent Dalitz plot analysis 

This analysis allows the measurement of a without am- 
biguities even in the presence of non-negligible penguin 
contributions. The B° — * amplitude is written as 

sum over the resonant intermediate states p + ft~, p~K + 
and p°n°: 

A = f + ~A+- +f- + A- + +/ 00 A 00 . (10) 

The form factors f + ~, f h , and / 00 describe the full kine- 
matics of the three-pion final state. They include the dom- 
inant p(770) as well as higher-mass states. An estimate of 
the form factors can be obtained from e + e~ — ► iz + n~ data 
and from T + — * 7T + 7r°V T decays. Toy Monte Carlo studies 
with samples of 500 events of pure signal (this corresponds 
to approximately 100 fb _1 ) with realistic At resolution and 
flavour tagging give typical statistical uncertainties on a of 
the order of 6° — 10° [ 16 1, depending on the values of the 
different amplitudes in Eq. (1101 . and for values consistent 
with the quasi-two-body result from BAB^R. In these toy 
fits, all (complex) amplitudes (T + ~, T h , P h and 
r 00 ) were free to vary and determined from the toy data. 
The inclusion of a realistic number of continuum back- 
ground events increases the statistical error by about 25 %. 
As in the quasi-two-body analysis, the impact of B-related 
backgrounds is important. Some of these backgrounds tend 
to accumulate in the interference regions of the Dalitz plot 
which are essential for the Dalitz plot analysis. The inclu- 
sion of realistic B-related backgrounds in the toys increases 
the statistical error by about 30 %. Shifts of up to 40° in 
the fitted value of a are observed in toy experiments where 
realistic B-backgrounds are simulated but not taken into ac- 
count in the fit 1 1161 . While B-backgrounds will be mod- 
elled in any real Dalitz plot analysis, the systematic uncer- 
tainties due to imperfections in the B-background model 
are potentially large. Various other issues will need to be 
addressed. These include the precise shape of the form fac- 
tors (moduli and phases). The impact of possible contribu- 
tions to B° — > 7i + 7i^7i° from very broad resonances and of 
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Figure 2. Expected constraints on a using SU(2) flavour symmetry for different scenarios (see text). 



Mode 



Observables 



2+0 A+0 

20+ 4 0+ 

200 ,^00 cOO 
5 p7r> ^pn' °pn 



^observables Experimental status 



B^ 

B+ 
B+ 
B° 



p°n 



6 
2 
2 
3 



Measured rfTTIITSl. 
Not observed. Limit on ^p®: 
Only 3$ p \ measured [EI|5l[T5l. 
Not observed. Limits on 0® n : [|3]|5]. 



Table 1. Observable-counting in the SU(2) analysis. 



= (22.6 ± 1.0) 10- 6 

= ( 0.9 ±0.3) 10~ 6 

= ( 8.1 ±1.0) 10~ 6 

0}+ = (8.7 ±0.5) 10~ 6 , 



A£p = -0.18 ±0.04, 

C p% = 0.36 ±0.08, 

AC pn = 0.28 ±0.08, 

S p7I = 0.19±0.10, 

AS pn = 0.15 ±0.11 , 



= 0.56 ±0.35, 

= 0.82 ±0.50, 

A$ = 0.02 ±0.07, 

A^p = 0.36 ±0.05, 



Table 2. Extrapolation to an integrated luminosity of 500 fb 
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non-resonant contributions needs to be evaluated. It could 
be possible to include a corresponding term in the ampli- 
tude i ll Ob and fit for its normalisation in the data. A more 
detailed discussion of these issues can be found, e.g., in [ 
E3[TU[T9l|20||2l]. This list of issues is not intended to be 
exhaustive. 

4 Conclusion 

The first time-dependent quasi-two-body analyses of B — > 
pn are becoming available from the B-factories. With 
SU(2) flavour symmetry, these measurements can be used 
to constrain the angle a. These constraints on a are only 
expected to be significant if the branching fraction of the 
yet unobserved decay B° — > p°7T turns out to be well be- 
low the current experimental sensitivity. The quasi-two- 
body analyses are also an important first step towards a 
time-dependent Dalitz plot analysis. This analysis could 
provide strong constraints on a with the luminosities ex- 
pected from the B-factories, but various experimental and 
theoretical issues still need to be addressed. 
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